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ABSTRACT 

GD 362 is a massive white dwarf with a spectrum suggesting a H-rich atmosphere which also shows 
very high abundances of Ca, Mg, Fe and other metals. However, for pure H-atmospheres the diffusion 
timescales are so short that very extreme assumptions have to be made to account for the observed 
abundances of metals. The most favored hypothesis is that the metals are accreted from either a 
dusty disk or from an asteroid belt. Here we propose that the envelope of GD 362 is dominated 
by He, which at these effective temperatures is almost completely invisible in the spectrum. This 
assumption strongly alleviates the problem, since the diffusion timescales are much larger for He- 
dominated atmospheres. We also propose that the He-dominated atmosphere of GD 362 is likely to 
be the result of the merger of a binary white dwarf, a very rare event in our Galaxy, since the expected 
galactic rate is ~ 10~ 2 yr -1 . 

Subject headings: stars: white dwarfs — stars: chemically peculiar — stars: individual (GD 362) 
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l2007h has shown that planets around white dwarfs with 
masses M\vd > 0.7 M are generally expected to be 
found at orbital radii r > 15 AU because they cannot 
survive the planetary nebula phase and that if planets 
are to be found at smaller orbital radii around mas- 
sive white dwarfs, they had to form as the result of 
the merger of two white dwarfs. It is also interesting 
to note that there have been previous suggestions about 
white dwarfs that are merger produc ts — see for instance 
iLiebert. Bergeron fe Holberg I (|2005f ) — but these claims 
do not have yet any observational support. 

2. THE SCENARIO 

Another possibility is that some massive white dwarfs 
are the result of the merger of a double white dwarf 
close binary system. This scenario has been stud- 
ied in several papers. However, in most of these pa- 
pers eit her the resulting nucleosynthesis was not ad- 
dressed (jSegretain. Chabrier. fe M ochkovitch 1993), or 
the spatial resolution was poor (|Benz et al. 111990 ). or 
the calcu lations were performed usin g crude approxi- 
mations (jMochkovitch fe Liviol Il990f ). Very recently, 
and using a Smoothed Particle Hydrodynamics code, 
a series of simulations with the adequate spatial reso- 
lution were performe d and the nucleosynthesis of the 
merger was studied dGuerrero. Isern. fe Garcia-Berro I 
l2004t lLoren-Aguilar et al. II2005D . The main results of 
such simulations are that the less massive white dwarf 
of the binary system is totally disrupted in a few or- 
„, , , , i t — 3 lu p. t jbital periods. A fraction of the secondary is directly 
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from the system. Ihe resulting temperatures are rather 



1. INTRODUCTION 

GD 362 has been interpreted as a massive, rather cool 
(T c ff k, 9740 ± 50 K), white dwarf with a heavy accretion 
disk surrounding it (iKilic et al~1 120051: iBecklin et al. I 
120051: iGianinas. Dufour. fe Bergeron I I2004D . The dusty 
disk around GD 362 produces an excess of infrared ra- 
diation which amounts to ~ 3% of the total stellar lu- 
minosity. The chemical composition of GD 362 is also 
rather singular, showing a hydrogen rich atmosphere 
wit h very high abundances of Ca, Mg, Fe a nd other met- 
als ((G ianinas. Dufour. fe Bergeron! 12004). Thus, it is 
classified as a massive DAZ (hydrogen-rich) white dwarf. 
The origin of such particularly high photospheric abun- 
dances — log(AW^H) = -5.2, log(Ar Mg /Ar H ) = -4.8 
and log^Fe/TVn) = —4.5 — and of the surrounding 
dusty disk around it still remains a mystery. Since the 
diffusion timescales for metals in H-rich white d warfs 
are of only a few years (|Koester fe Wilkeril [20061) very 
extreme assumptions have to be made in order to ex- 
plain these abundances. At present the most widely ac- 
cepted scenario is disruption and accretion of a plane- 
tary body, although for this scenario to be feasible the 
planetary system should survive during the advanced 
stages of stellar evolution, which by no means is guaran- 
teed. Thus, the formation of an asteroid would require 
the previous existence of a disk around this white dwarf 
(ILivio. Pringle. fe Saffer 1 1 19921 iLivio. Pringle. fe Wood] 
2005). Particularly, a recent analysis (|Villaver fe Liviol 



2 



Garcfa-Berro et al. 



TABLE 1 

Main results of the SPH simulations. 



Run 0.4+0.8 0.4+1.0 0.4+1.2 0.6+0.6 0.6+0.8 



Af WD /Mq 099 1+6 UO" O90 +09" 

M disk /M e 0.21 0.24 0.30 0.30 0.29 

A/ ci /M Q 10~ 3 10~ 3 10~ 3 10~ 2 10~ 3 

He 0.94 0.93 0.99 

C 3 x 10~ 2 2 x 10~ 2 5 x 10~ 3 0.4 0.4 

O 1 x 10~ 2 3 x 10~ 3 3 x 10~ 3 0.6 0.6 

Ca 4 x 10~ 5 2 x 10~ 4 9 x 10~ 6 

Mg 3 x 10~ 5 3 x 10~ 5 6 x 10" 6 

S 8 x 10 -5 2 x 10" 4 5 x 10" 7 

Si 1 x 10~ 4 2 x 10" 4 3 x 10" 5 

Fe 9 x 10~ 3 7 x 10" 3 5 x 10~ 4 



high (~9x 10 s K) during the most violent phases of the 
merger, allowing for extensive nuclear processing. 

The enhancement of the abundances of the most rele- 
vant nuclear isotopes occurs when one of the coalescing 
white dwarfs is made of pure He. Table 1 shows the 
average chemical composition of the resulting disk and 
the main characteristics of some selected simulations. It 
should be noted, however, that the distribution of the 
different elements in the disk is rather inhomogeneous. 
Obviously those parts of the disk in which the mate- 
rial of the secondary has been shocked have undergone 
major nuclear processing. Hence, these regions are C- 
and O-depleted and Si- and Fe-enhanced. In fact, the 
innermost regions (R < 0.1 Rq) of the merged object, 
which have approximately the shape of an ellipsoid, are 
C- and O-rich. It is expected that this region would be 
eventually accreted during the the first moments of the 
cooling phase of the central object, leading to a more 
massive white dwarf. We also find that the abundance 
of intermediate-mass and iron-group elements is consid- 
erably larger than th at of C and O in the remnants of 
the a ccretion stream (jGuerrero. Isern. fe Garcfa-Berro 1 
l200l which are at larger distances, thus favoring smaller 
accretion rates in order to explain the Ca abundance. In 
any case, if the photospheric abundances of GD 362 are 
to be explained with this scenario the accretion of He 
rich material is required. 

Since He is also accreted onto the surface of GD 362, 
the photospheric layers may contain significant amounts 
of He which, at the effective temperature of GD 362 
would be almost spectroscopically invisible. Thus, 
GD 362 would still be classified as a DA white dwarf 
provided that some H is present in its atmosphere. Con- 
sequently, the H/He ratio can be regarded as a free 
parameter. However, the presence of He in a cool 
hydrogen-rich atmosphere affects the surface gravity de- 
termined from spectroscopy, and thus the mass d eter- 
mination (Ber geron. Wesemael. fe Fontainel |1991[ ) . In 
Fig. 1 we show three almost identical synthetic spec- 
tra representative of GD 362 with various assumed He 
abundances. If He/H=10 is adopted then logg = 8.25 is 
obtained (Mwd ~ 0.8 Mq) whereas if we adopt He/H=l 
then the surface gravity turns out to be log<? = 8.72. 
This corresponds to a mass of the primary of Mwd ~ 
1.0 Mq, which can be obtained from the coalescence of 
a 0.4 + 0.8 Mq binary system. Additionally, in this case 
the largest abundances of the relevant elements are ob- 
tained. Thus, we choose the 0.4 + 0.8 Mq simulation as 
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Fig. 1. — Spectrum of GD 362 for three different helium 
abundances. The black line shows the spectrum of GD 362 
when a pure hydrogen atmosphere is assumed, leading to a sur- 
face gravity of logg = 9.12. For increasing amounts of He — 
namely A+He)/A+H)=l, red curve, and A r (He)/A r (H)=10, blue 
curve — the corresponding surface gravities are smaller. The inset 
shows an expanded view of the predicted He line at 5876 A for 
A r (He)/A r (H)=10. High quality spectroscopic observations should 
be able to confirm its presence, which has been recently reported 
llJura et al. 11200^) . See the electronic edition of the Journal for a 
color version of this figure. 

our reference model, although reasonable results can be 
obtained adopting other masses. In passing, we note that 
nevertheless the He abundance is rather uncertain since 
equally good fits to the observed spectrum of GD 362 can 
be obtained with very different He abundances. Thus, 
the mass of GD 362 is also rather uncertain. More im- 
portantly, if the mechanism producing the unusual pho- 
tospheric abundance pattern of GD 362 were to be ac- 
cretion from the inner regions of the disk — which are 
C- and O-rich — atomic lines of CI, and C2 molecular 
bands should be rather apparent in the spectrum. But 
the strength of these spectral features depends very much 
on the adopted He abundance, because the atmospheric 
pressure and the opacity also depend very much on the 
H/He ratio, which is rather uncertain. 

In order to know whether the chemical abundances 
of GD 362 can be reproduced by direct accretion from 
the keplerian disk we proceed as follows. Given the sur- 
face gravity and the effective temperature of our model 
we compute the luminosity, the radius and the cooling 
time of th e white dwarf acco rding to a set of cooling 
sequences (IS alar is et al. I [2000h . We obtain respectively 
log(L WD /L ) -3.283, \og(R WD /R Q ) ~ -2.096, and 
t coo \ ~ 2.2 Gyr. Hence, in this scenario GD 362 has 
had enough time from the moment in which the merger 
occurred to cool down, to accrete most of the C- and 
O-rich region, settle down the accretion disk, and to 
form dust. Additionally, the central white dwarf has had 
time enough to accrete (at a rate much smaller than the 
Bondi-Hoyle accretion rate) the small amount of hydro- 
gen from the ISM to show spectroscopic hydrogen fea- 
tures. We further assume that the accretion luminosity: 

i» - (i) 

HwD 

is, in the worst of the cases, smaller than the luminosity 
of the white dwarf. This provides us with an (extreme) 
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Fig. 2. — Spectral energy distribution of GD 362. The figure 
shows the spectral energy distribution of GD 362. The dotted line 
shows the spectrum of a white dwarf with and effective temperature 
of 9740 K and logg = 8.72, which corresponds to a mass of about 
one solar mass, the dashed line shows the spectrum of a passive fiat, 
opaque dust disk and the solid line depicts t he composite sp e ctrum . 
The observational data were obtained from Bccklin ct al. (2005). 



upper limit to the accretion rate, which turns out to be 
1.3 x 10 -13 M yr -1 . Next, we assume that the abun- 
dance of Ca is the result of the equilibrium between the 
accreted material and gravitational diffusion: 



Tdiff 



(2) 



where -X"di s k is the abundance in the accretion disk, 
X Q bs is the photospheric abundance, M cnv is the mass 
of the envelope of GD 362 and Td;s is the diffusion 
timescale. The diffusion timescale of Ca for H-rich at- 
mospheres is of the order of a few years. However, the 
accreted material is He-rich, so the diffusion timescale 
is probably mo re typical of a He-rich envelope, which 
is much larger (|Paauette et al. |[l986h . of the order of 
Tdiff ~ 1.5 X 10 4 yr. Unfortunately, diffusion timescales 
for mixed H/He envelopes do not exist. However, the 
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diffusion characteristic tim es scale as Tdiff °c pT~ 
(|Alcock fc Illarionov ll!980T ). We have computed detailed 
atmosphere models for pure H, He/H=l and He/H=10 
and scaled the diffusion timescale using the values of 
the density and the temperature at the base of the 
convective zones and the appropriate chemical composi- 
tion. For our fiducial composition (He/H=10) we obtain 
Tdiff ~ 8.5 x 10 3 yr. From this we obtain the mass of 
the region where diffusion occurs, which turns out to be 
M env ~ 7.2 x 1CP 9 Mq, which is much smaller than that 
obtained by accretion from the interstellar medium at the 
Bondi-Hoyle accretion rate (~ 1.5 x lCP 6 Af Q ). Hence, 
the photospheric abundances of GD 362 can be success- 
fully explained by direct accretion from the surrounding 
disk. 

Now we assess whether the flux from the accretion disk 
can be fitted by the results of our SPH simulations. In 
order to compute the flux radiated away from the system 
two contributions must be taken into account. The first 
one is the expected photospheric flux from the star, for 



Fig. 3. — Evolution of the rotational velocity for several field 
strengths, the observational upper limit is shown as a horizontal 
dashed line. 



which we use the spectral energy distribution (Bwd) of 
a white dwarf of mass 1 Mq, at T e g w 9740 K: 



i?wr 



WD 



-BwD(Icff) 



(3) 



Given the luminosity of our model and the apparent 
magnitude of GD 362 we obtain a distance of Z>wd = 
33 pc. The second contribution to the total flux comes 
from the emission of the disk, which for a passive flat, 
opaq ue dust disk is (|Chiang fc Goldreich I Il997t Uura I 
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where i is the inclination of the disk (which we adopt to 
be face-on), x m — hvjk^T-^ and T- m = 1200 K is the 
condensation temperature of silicate dust. The outer ra- 
dius is taken from the results of our SPH simulations and 
turns out to be i? ut ~ 1 Rq- The result is displayed in 
Fig. 2. The dots are the observational data for GD 362. 

The proposed scenario has apparently two weak points. 
The first one is that infrared observations indicate the 
presence of SiO. This requires that O should be more 
abundant than C in order to form it. However our simu- 
lations show that the ratio of C to O is a function of the 
distance to the primary and, in some regions of the disk 
the ratio is smaller than 1 , allowing for the formation of 
SiO in the accretion disk. Furthermore, after 2.2 Gyr of 
evolution the resulting disk has had time to form planets 
or asteroids with the subsequent chemical differentiation. 

The second apparent drawback of the model is that 
the central white dwarf rotates very fast. However, an 
unobservable magnetic field can brake down the central 
star to acceptable velocities. Using the observed spec- 
trum of GD 362 it is possible to set an upper limit to 
the rotation speed of usini < 500 km s _1 . We assume 
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that the central white dwarf has a weak magnetic field, 
B. The magnetic torques that lead to spin-down are 
caused by the interaction between the white dwarf and 
the surrounding disk. The evolution of the angular ve- 
locity due to the coupling of the white dwarf magneto- 
sphere and the disk is g i yen bv (iLivio fc Pringle 1 1991 
lArmitage fe Clark I [1999 : iBenacouista et al. 112003 ): 



_V^3 

~~ He? 



2 

sm a 



El. 

3/ 



1 



(5) 

where \i = i3i?| VD , R m is the magnetospheric radius of 
the star, / is the moment of inertia, a is the angle be- 
tween the rotation and magnetic axes (which we adopt 
to be 30°) and 



R c 



I GA/wd 



1/3 



(6) 



is the corotation radius. The first term in this expression 
corresponds to the magnetic dipole radiation emission, 
the second to the disk-field coupling and the last one 
to the angular momentum transferred from the disk to 
the white dwarf. The magnetic linkage between the star 
and the disk leads to a spin-down torque on the star if 
the magnetospheric radius is large enough relative to the 
corotation radius: 



| > 2 _ 2 /3 

R c 



(7) 



We adopt R m = R c . Solving numerically the previous 
differential equation with the appropriate parameters for 



our case, the evolution of the rotation velocity is shown 
in figure 3. As can be seen a weak magnetic field of 
about 50 kG is able to brake down the white dwarf to 
velocities below the observational limit. This magnetic 
field is much smaller than the upper limit of about 0.7 
MG obtained from the spectrum of GD 362. Hence, our 
scenario also accounts for the low rotational velocity of 
GD 362, without adopting extreme assumptions. 

3. CONCLUSIONS 

We have shown that the anomalous photospheric chem- 
ical composition of the DAZ white dwarf GD 362 and of 
the infrared excess of surrounding disk can be quite nat- 
urally explained assuming that this white dwarf is the 
result of the coalescence of a binary white dwarf system. 
This scenario provides a natural explanation of both the 
observed photospheric abundances of GD 362 and of its 
infrared excess without the need to invoke extreme as- 
sumptions, like the accretion of a planet or an asteroid, 
since direct accretion from the disk surrounding disk pro- 
vides a self-consistent way of polluting the envelope of 
the white dwarf with the required amounts of Ca, Mg, 
Si and Fe. Moreover, this last scenario can be also well 
accomodated within the framework of our scenario given 
that the formation of planets and other minor bodies is 
strongly enhanced in metal-rich disks. Hence, GD 362 
could be the relic of a very rare event in our Galaxy: the 
coalescence of a double white dwarf binary system. 
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